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I. INTRODUCTION 


In this report the geometrical theory of diffraction (GTD) [1] is 
used to treat the diffraction by a perfectly-conducting rectangular 
cylinder which is illuminated by an array of electric or magnetic current 
line sources. The illumination by a single line source of the array is 
considered first with the geometry of the problem shown in Fig. 1. The 
scattered and total fields are calculated at all points exterior to the 
rectangular cylinder, except for the shaded regions around the source 
and edges, which are excluded because of the nature of the high-frequency 
approximation. More will be said about this later. The solution is 
then generalized to calculate the field of linear arrays of such line, 
sources radiating in the presence of the cylinder. 

The scattering by a rectangular cylinder has been considered pre- 
viously for the special case of plane wave incidence. Mei and VanBladel 
1 2 , 3 J formulated the problem in terms of an integral equation which they 
solved numerically to find the surface currents on the rectangular 
cylinder. They give radiation patterns, scattering cross sections, and 
surface currents for both E- and H-polarized incident waves. There are 
two disadvantages associated with this method: 1) as the frequency 

increases, difficulties are encountered with the convergence of the 
solution (this is particularly true when the technique is employed to 
solve a three-dimensional scattering problem), 2) the solution provides 
no physical insight into the scattering mechanism. Morse |4] also 
studied this problem using the ordinary GTD to obtain expressions for 
the diffracted field away from shadow and reflection boundaries. Since 
the ordinary theory fails at shadow and reflection boundaries, he intro- 
duced supplementary solutions there. He employed Oberhettinger' s uni- 
form asymptotic solution J 5 ] near the boundaries of the incident and 
reflected fields, and he employed an integral representation of the 
field near the shadow boundaries of the fields of the diffracted rays. 
Thus, he did not obtain a compact high-frequency solution to this problem 
The difficulties encountered by Morse at these boundaries may be overcome 
with a new edge diffraction coefficient derived by Kouyoumjian and Pathak 
|6,7|. This diffraction coefficient can be applied in the transition 
regions adjacent to the shadow and reflection boundaries so that one ob- 
tains a total field which is valid and continuous everywhere away from 
edges and caustics. Their diffraction coefficient is employed in this 
analysis . 

Using Keller's Generalized Fermat's Principle, we include contri- 
butions to the total field from the geometrical optics fields (incident 
and reflected), as well as singly-di ffracted fields which appear to 
emanate from the edges. Doubly-diffracted fields have also been in- 
cluded to describe the interactions between the edges; however, multiply- 
diffracted rays of higher order have been neglected because, in general, 
their fields contribute insignificantly, for the problem defined here. 

The incident, reflected and diffracted rays may be shadowed in the geo- 
metrical optics sense, and hence they contribute to the total field only 
in their respective regions of illumination. Discontinuities in the 




Fig. 1. Line source in the presence of the rectangular 
cylinder. 


field are introduced at shadow boundaries and at the reflection boundaries 
but they are symstematically removed by employing the new diffraction co- 
efficient for the edge diffracted field mentioned previously. It should 
be pointed out that there are some residual discontinuities due to uncom- 
pensated discontinuities in the field of the doubly-diffracted rays at 
their shadow boundaries. However, these discontinuities are so small 
that they are not apparent in the plotted patterns. The patterns calcu- 
lated from our GTD solution are found to be in excellent agreement with 
those calculated from numerical solutions. 

Due to the high-frequency approximations of the incident and dif- 
fracted fields, our solution is restricted so that the following distances 
are greater than 0.7 wavelength. 

1) The distance between the line source and the closest edge 
of the rectangular cylinder, 

2) the distance between the observation point P and the 
closest edge of the rectangular cylinder, 

3) the distances between the edges of the rectangular 
cylinder, 

4) the distance between the observation point P and the-, 
closest line source, when calculating the incident. or 
total fields. 
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A computer program based on the GTD solution for the rectangular 
cylinder in the presence of an array source has been written and is 
included in the report. With this program, it is possible to obtain 
numerical results for both the near- and far- fields of the cylinder under 
quite general conditions of illumination. Thus, the program is directly 
relevant to both antenna and scattering problems. As an example of its 
versatility, one notes that the program may be used to compute the 
pattern of an array of magnetic line sources mounted directly on the 
rectangular cylinder,, provided the sources are not too close to an edge. 

II. METHOD OF SOLUTION 

Keller's geometrical theory of diffraction j'l] is an extension of 
geometrical optics in which diffracted rays are introduced by a generali- 
zation of Fermat's principle, the excitation of the diffracted field is 
treated as a local phenomenon, and away from the diffracting surface 
the behavior of the diffracted field along its ray is the same as that 
of the geometrical optics field. The basic idea of GTD is that the 
field of the line source illuminates the rectangular cylinder giving 
rise to a reflected field and an edge diffracted field, which consists of 
the fields of singly and multiply-di ffracted rays. The total field U(P) 
at a point P is equal to the sum of the fields on all rays through P. 


(1) U(P) = I U-j ( P ) 

rays 

which includes the incident field if P is not in the shadow region. The 
wave function U(P) represents a magnetic field parallel to the edge in 
case of a magnetic line source, and an electric field parallel to the 
edge in the case of an electric current line source. The pertinent rays 
and their associated fields will be discussed briefly in the following 
paragraphs . 

A. Geometrical Optics Rays - Direct and Reflected 

Let us consider the field radiated from a line source at 0 and 
observed at P as shown in Fig. 2. Fermat's principle predicts only the 
direct ray OP. If a line source is being considered, the field along OP 
is given by 

-jks 0 

(2) U 1 (P) = C - — — 

/s 0 

where s Q is the distance between 0 and P, and C is a conveniently chosen 
normalization constant. For the configuration shown in Fig. 2, the space 
surrounding the right-angle wedge may be divided into three regions: 
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SHADOW BOUNDARY AT tt '+ <j>' 
REFLECTION BOUNDARY AT tt - <f> 


Fig. 2. Line source in the presence of a wedge. 


Region I 0 < < 'tt - <j> 1 , 

Region II tt - ^> ’ < '4> < % + <(.', 

Region III ir + <t>‘ < <p< 3ir/2. 

Region III is the shadow region, which is not penetrated by the incident 
ray; the incident field vanishes here. 

We know that there is a field reflected from the surface AQr. To 
describe this we introduce an additional class of rays which include on 
their trajectory a point Qr of the surface AQ^. Applying Fermat's 
principle, the distance OQrP along the ray path is a minimum and the law 
of reflection results. This simple extension of Fermat's principle which 
accounts for the reflected ray is so natural that we accept it without 
question. The field of the reflected ray is readily deduced from image 
theory as 


-jks" 

(3) U r (P) e U r (A-Q E ) = ± C , 

where the positive sign is for Neumann (hard) boundary condition associated 
with the magnetic current line source, the negative sign is for the 
Dirichlet (soft) boundary condition associated with the electric current 
source line, and s" is the distance between the image O' and the observa- 
tion point P. The reflected field vanishes in regions II and III, which 
the reflected ray does not penetrate. Let us consider now a further ex- 
tension of Format's principle. 
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B. Si n gl y-Di ff racted Rays 

It is well known that the ray incident on the edge Q E in Fig. 2 
gives rise to diffraction. To account for this, Keller introduced a 
class of rays which includes the point Qr in its trajectory. This com- 
pletely determines the diffracted ray path in the isotropic, homogeneous 
medium of this two-dimensional problem, so the law of edge diffraction 
becomes trivial under these circumstances. 

In terms of GTD, the diffracted field at P for the line source at 

0 is 

~j ks 

(4) U d (Q E ) = U^Qe) Dt(M') ~z — , 

n 


where is the scalar diffraction coefficient for the accoustically soft 
(Dirichlet) boundary condition and D h is the scalar diffraction coeffi- 
cient for the accoustically hard (Neumann) boundary condition. They are 
deduced from the general dyadic diffraction coefficient D((f>,<f>' ,e 0 ) ob- 
tained by Kouyoumjian and Pathak [6,7]. For the special case where the 
incident ray is perpendicular to a straight edge, the scalar diffraction 
coefficients are given by 


(5) 


Ds ( <f>» ) 

h 


_ e -J */4 

2n^2irk 


f fn + (<P-<P') 

COtl 


J F[kLa + (4> -4> ' )J + cotj — — J F [kL a" (<)>-<}>' )] 


2n 


cot 


nr+ (<(> +< f> 1 ) 


)\ /TT- ( 4 >+ci> ' )\ 

jF[kLa (<j>+<j>')] + cotf 2n ) F[kLa" ( cj>+ 4 > ' )] 


2n 


where ntr is the exterior wedge angle, which equals 3ir/2 in this case, and 


(6) F(x) = 2d|^3r| e jX J e" jT dr 

I.AI 

in whi ch 


(7) 


a ± ( 4>±4> ' ) 


fcnirN 1 
2 cos 2 — — 


2 / ' 
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N + are the integers which most nearly satisfy the following equations 

(8) 2irnN = 7r + (<f>±4>') 

(9) 2-n-nN = — ft + ( 4>±4> 1 ) 


and kL is the large parameter in the asymptotic evaluation of the pertinent 
integrals involved in the derivation of the dyadic diffraction coeffi- 
cient. The quantity L may be viewed as a distance parameter whi ch depends 
upon the type of edge illumination; for line source illumination, L is 
given by 


s s 1 

do) i ■ S -TT 

For grazing incidence (f>' = 0, ott, Dh is multiplied by a factor of 
furthermore, if the diffracted ray grazes the surface in the case of 
a soft boundary, D = 0 and the diffracted field vanishes, as it 
should. 

the field of the singly-diffracted ray is discontinuous at the 
shadow and reflection boundaries in a way which compensates the dis- 
continuities in the geometrical optics fields there. This is readily 
demonstrated; consider' for example the incident and diffracted fields at 
the shadow boundary, where to simplify the discussion, it is assumed 
there is no nearby reflection boundary. Let tt+^'-e be a point close to 
the shadow boundary, see Fig. 2. In the illuminated region e > 0 and in 
the shadow region e < 0. 


(ID 


UU) =< 


C 


-j' ks o 

e 




+ C 


-jks ' 
e 




-jks 

e 

D S ( 4> ' +TT-E , 4> 1 ) 

h viT 



D^U'+TT-ejtjj' ) 


e 


-jks 


Js 


E < 0 


E > 0 


For e small it follows from Eq. (5), that 


-e 


-j tt/4 


(12) D s (4>'+tt-e,4) 1 ) - 

n 2nv#7rk 


cot ^ F[kLa"U-E)] 


+ smaller terms which are continuous at the shadow boundary | . 
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From Eq. (9), 


N" = 0 
Also, as e -v 0, 

(13) cot (§-) = — 

2n e 

( 14 ) a'U-e) = ~ 

jir/4 

(15) F[kLa‘(ir-e)] = /d<L72 e |e| 

Substituting Eqs. (13) and (15) into Eq. (12) as e -> 0, 

i 5 1 s 

(16) D s (<t> , +TT-e,<j>' ) = -h\ ; sgn e + smaller, continuous terms, 

h v s + s 

(17) s 0 = s 1 + s 


Upon substituting Eqs. (16) and (17) into Eq. (11), it is seen that 
the total field is continuous at the shadow boundary. In an analogous 
manner it can be shown that the total field is continuous at the reflec- 
tion boundary. 

C. Doubly-Diffracted Rays 

When one face of the conducting wedge is terminated at Qp as shown 
in Fig. 3, a second order di ffracted-ray will emanate from the edge Qp. 
In terms of the GTD, the doubly-diffracted field at P due to the line 
source at 0 can be written as 

“ jks 

(18) U d (Q E ,Q F ) = UMQ f ) h Dc ( 4> 2 >0 ) 

n i/s 


3tt 


= ju 1 (q e ) d^(|- ,♦') 


j-jkh| 


\ 8c ( 4 ) 2 >8 ) 


S -Jks 

T" 
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Since D s ( 3 tt/2 , 4 > 1 ) = 0, the contribution from the doubly-diffracted rays 
vanishes for the soft boundary according to the above expression. If a 
higher order approximation for the doubly-di ffracted field is employed, 
then this contribution is non-vanishing, as will be explained later. 

The field of the ray singly-diffracted at has a shadow boundary 
SB(Qp) as shown in Fig. 3; the singly-di ffracted ray does not penetrate 

0 



the shaded region. It will be shown next that the discontinuity in the 
field of the singly-diffracted ray at SB(Qp) is compensated by the ray 
doubly-diffracted from Qp, so that the total diffracted field is con- 
tinuous at this boundary. Since the field doubly diffracted vanishes 
in the case of the soft boundary, we only need to treat the hard 
boundary here. 


Consider a point close to SB(Q E ) so that <j> 2 = n-e, the total 
diffracted field at this boundary is 


(19) U TD = { 


u 1 (Q e ) yv*') uMQp) 


-J' ks l Mir-e.O) -jks 


Js 


, e > 0 


. D h (w-e,0) -jks 

( U (Qc) -5-r , e< 0 

c Js 
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where s 1 ,<j> 1 are the coordinates of the ray diffracted from Q E , and 

(20) U 1 (Q F ) = U 1 (Q E )D h (|2- ,*'} . 

When <j >2 = 0 and the singly-diffracted ray grazes the vertical surface, 
the second and fourth terms in the expression for the diffraction co- 
efficient are the same, except for the ± sign of the latter. This is 
also true for the first and third terms. As a result. 


_ p -j tt/4 

(21) D h (ir-e,0) = — 

nv^Trk 

+ smaller terms which are continuous at SB(Q E ) 


As e -> 0, it is seen from Eqs . (12), (13), (14) and (15) that 


(22) D h (u-e,0) = - 


furthermore. 


hs 

1+h s 9 n e ’ 



jcot(|^-) F[kLa"(7T-e)] 


(23) s x = h + s 


Substituting Eqs. (22) and (23) into Eq. (19) and making use of 
Eq. (20), it is seen that the total diffracted field is continuous at 
the boundary SB(Qjr). 


As we have already noted, in the case of a soft boundary the field 
of an incident-ray grazing the surface vanishes, the edge-di ffracted 
field is then proportional to the normal derivative of the incident 
field at the edge. The proportionality factor is a diffraction coef- 
ficient D 1 given by Karp and Keller |8J. Thus, for the case of Dirichlet 
problem, the doubly-diffracted field must be replaced by 


(24) 


d aUMQr) D 1 (<j> 2 ,0) -jks 

U(Q E ,Qr) = — £ 

h F an 2 \/s 


where 



The derivate aU(Qp) / an is taken with respect to the normal to the 
surface QfQf* This contribution is weak in comparison with that of the 
singly-diffracted rays; the. contribution of the former is of order (1/k 2 ), 
whereas that of the latter is of order (l/vIT). In calculating the field 
diffracted from the soft cylinder, it was found that the field of the 
doubly-diffracted rays did not contribute significantly, so the contribu- 
tion from these rays can be omitted in this case. 

Let us now turn to the diffraction by a rectangular cylinder 
illuminated by a line source. Depending on the location of the line 
source, the whole domain surrounded by the cylinder will be divided into 
regions by the various shadow boundaries and the reflection boundaries. 
Each of these boundaries is labeled to indicate how it originates. For 
example, referring to Fig. 4, the notation^SB means the shadow boundary 
of the incident geometrical optics field U' r (P), RB(A-B) is the shadow 
boundary of the geometrical optics field U r (A-B) reflected from the . 
surface A-B, SB (A) is the shadow boundary of the singly-diffracted U d (A), 
which emanates from the edge A, and SB(A,B) is the shadow boundary of 
the doubly-diffracted field U a (A,B) which emanates from the edge B. The 
shadow boundary of the reflected field is referred to simply as the 
reflection boundary. 


R.B.( A-B.) S.B.(A,D.) 



S.B.( A.) 
S.B.(B,A.) 


S.B. ( B. ) 
S.B.(A,B) 


Fig. 4. Shadow and reflection boundaries of the 
GTD fields. 
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CHECK (/) MEANS THE FIELD IS NON- VANISHING IN THE REGION 
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In what follows we may locate the line source in the upper left 
hand quadrant of the space surrounding the cylinder, without loss of 
generality. Three cases will be considered. 

1) First, the line source is located directly above the cylinder. 

The whole domain around the cylinder is divided into 14 regions 
as shown in Fig. 4. Table I shows the regions covered by each 
type of ray. To demonstrate the use of the table, let us 
consider an observation point P in region IX. Examining the 
column under region IX of Table I, one finds checks (/) for 
U d (A), U d (B,A), U d (A,D), and U d (B,C). Thus the total field 

at P is equal to U T (P) = U d (A) + U d (B,A) U d (A,D) + U d (B,C). 

2) The line source is located to the upper left of the cylinder 
as shown in Fig. 5. In this case there are 18 regions. 

Table II shows the regions where a particular ray exists and 
its field makes a non-vanishing contribution. 

3) Line source at grazing incidence. Let the line source be 

located in the plane which contains the A-D side of the rec- 
tangular cylinder as shown in Fig. 6. This case is of special 

interest because the trailing edge D lies in the shadow of 

the leading edge A. Also, it lies on the shadow boundary of 

the direct geometrical optics field. The behavior of the 
field at the boundary DP must be treated separately, see 
Appendix I. As before, the domain surrounding the cylinder is 
divided into regions as shown in Fig. 6 and Table III gives the 
regions covered by the individual rays. 

Recall that one or more of the various field components is dis- 
continuous at each boundary shown in Figs. 4, 5 and 6 but that all except 
the discontinuities in the doubly diffracted fields are compensated; e.g., 
the discontinuities in the geometrical optics field are compensated by 
the field of the singly-diffracted rays and the discontinuities in the 
field of the singly-diffracted rays are compensated by the fields of the 
doubly-diffracted rays. 

As a final step in the anslysis, the fields of the individual line 
sources are superimposed to give the field of a linear array of line 
sources radiating in the presence of the cylinder, see Fig. 7. A com- 
puter program has been written to calculate the incident, total and 
scattered fields once the linear array is specified. Unlike its earlier 
definition for geometrical optics, the term incident field used here 
means the field of the array in the absence of the cylinder, and the 
scattered field is simply the difference between the total field and 
this incident field. The versatility of such a program is evident; the 
scattering from the cylinder for a wide variety of illuminations can be 
studied, and the radiation from antennas in the presence of the rectan- 
gular cylinder also can be studied. As a matter of fact, the program 
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was written originally so that the linear array of line sources, when 
densely packed, closely approximates the field of an aperture antenna of 
finite width W. The aperture antenna (more precisely its axis) is 
directed toward a point Q on the surface of the rectangular cylinder as 
shown in, Fig. 7. A description of the aperture radiation in terms of an 
array of discrete line sources is discussed in the following paragraphs. 

The aperture distribution may be approximated by a discrete array 
of line sources which are properly weighted in amplitude and phase. The 
width of the aperture denoted by W, is divided into 2M segments; (M = 
integer). The line sources are positioned at the ends of these segments, 
which introduces 2M + 1 line sources. In approximating a continuous 
distribution, the number M is selected so that 2M + 1 10 W/x, where 

X = free space radiated wavelength. 

Three types of line sources are available in this program: 

Type I An electric current line source 

Type II A magnetic current line source 

Type III A magnetic current moment line source. 

As described earlier, the electric current line source radiates an 
omnidirectional electric field which is parallel to the edge of the 
rectangular cylinder, and the magnetic current line source radiates an 
omnidirectional magnetic field which is parallel to the edge of the 
rectangular cylinder. The magnetic current moment line source consists 
of a continuous array of magnetic current moments directed perpendicular 
to the line of the array and parallel to the aperture in question. This 
line source radiates an electric field parallel to the edge of the 
cylinder; however, the field has a pattern, |.cos e|, where e is shown in 
Fig. 7. The strength of these magnetic type line., sources is determined 
irom^the equivalent magnetic surface currents in the aperture. K s = 

E x n, where E is the electric field distribution in the aperture 
(assumed known) and n is the outward normal to the aperture. 

The field of the two-dimensional aperture can be adequately repre- 
sented in the forward direction by a densely-packed array of type II and 
type III line sources, but such an array fails to approximate the field 
adequately at aspects behind the aperture. This limitation is particu- 
larly troublesome when calculating the total field. To overcome this 
difficulty an obliquity factor has been included in the program which 
multiplies the pattern of each line source. The obliquity factor is 
f(e) = cos n 0/2, where n = 0, %, 1, 2. When n = 0, the obliquity factor 
is unity so that the array radiates symmetrically with respect to the 
axis of its elements. The case n = 2 occurs naturally in the description 
of the radiation fields of aperture antennas via the Kirchhoff-Huygen's 
approximation (for the forward region ). The cases n = %, and n = 1 are 
added so that the n which best approximates the measured aperture pattern 
may be used. It is evident that the obliquity factor results in a 
pattern null in the direction directly behind the aperture at e =tt. In 
most practical cases, there is no such null in the backward direction. 
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Shadow and reflection boundaries of the 6TD fields 
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TABLE II 
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adow and reflection boundaries of the GTD fields. 
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SEE APPENDIX 
















































































2-0 APERTURE R i 
REPLACED BY AN 
ARRAY OF 2-D 
SOURCES 
( APERTURE 

WIDTH = W ) 


(OBSERVATION 
POINT ) 


Fig. 7. An array of line source in the presence of 
the rectangular cylinder. 


On the other hand, the pattern in the forward region is quite satisfac- 
tory when the obliquity factor is included, and as one moves away from 
the forward region, the pattern drops to a level where the differences 
between the simulated and measured patterns are unimportant when com- 
puting the total field surrounding the rectangular cylinder. 

The scattered, and the total fields are computed as a function of 
9s(-tt <_ ©s ir ) at a range R s . The incident field pattern on the other 
hand is available for a phase reference at the center of the array 
(aperture) at 0 in Fig. 7 and computed at a range measured from 0, as 
well as for a phase reference at Q with a range R s (just as for the 
scattered and the total fields). 

The input and output variables for the computer program are des- 
cribed in Appendix II, and the listing of the computer program is given. 
Also, a sample case is treated in this appendix to illustrate the use 
of the program. 

III. NUMERICAL RESULTS 

To assess the accuracy of the GTD solution described injthe 
preceding section, it was applied to several simple examples-, where the 
rectangular cylinder is illuminated by either an electric o!^ magnetic 
current line source, Far-zone patterns for the total field are calcu- 
— lated by this method and al^o from numerical solutions of the pertinent 
i n tegraU-equation *. The cylinders are square with a side length of 1.6 


*The computer programs for the integral equation solutions were provided 
by Prof. J.H. Richmond of the ElectroScience Laboratory. 
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wavelength, and in each case the line source illuminates the cylinder 
from a distance of 0.8 wavelength. These small distances provide a 
stringent test of our GTD solution; also they give us an opportunity to 
examine the accuracy of the new scalar diffraction coefficients in a 
situation where the edges are illuminated by curved wavefronts and where 
the transition regions are relatively broad. The pattern calculated 
from the integral equation solution must be considered more accurate 
for the small dimensions chosen for these examples, since the integral 
equation method is convergent whereas the GTD solution is an asymptotic 
approximation. 

Patterns for magnetic current line source excitation (hard boundary 
case) are given in Figs. 8, 9 and 10, where the line source is positioned 
on the diagonal of the square cylinder, on the centerline directly above 
the cylinder, and at a point of glancing incidence on one of its sur- 
faces. The agreement between the patterns calculated by the GTD and 
the integral equation method is remarkable - every detail is the same 
within the limits of graphical accuracy. The corresponding patterns for 
electric current line source excitation (soft boundary case) are present- 
ed in Figs. 11, 12 and 13. Again there is excellent agreement between 
the two pattern calculations, except in the vicinity region of forward 
scatter in Figs. 12 and 13. Note that the level of the patterns is very 
low in these regions, so that small errors in the solution become signi- 
ficant. We hope to look into the reason for these differences at a 
later time. 

' The numerical examples considered here confirm the accuracy and 
applicability of our GTD solution; this is further demonstrated by an 
example treated in Appendix II. As the size of the cylinder and the 
distance between the edges of the cylinder and the source (or sources) 
increases in terms of a wavelength, one can expect the accuracy of the 
GTD solution to increase, because it is an asymptotic approximation where 
k = 2irA is a large parameter. 

In the case of magnetic current line source excitation there is 
little evidence of shadowing in the forward direction by the small square 
cylinder; however, there is distinct evidence of shadowing in the case 
of electric current line source excitation, where the total electric 
field is parallel to the cylinder and must vanish at its surface. One 
should expect this . 
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6TD 



180 ° 


Fig. 8. Pattern of a magnetic current line source in 
the presence of a rectangular cylinder. 
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Fig. 9. Pattern of a magnetic current line source in 
the presence of a rectangular cylinder. 
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Fig. 11. Pattern of an electric current line source in 
the presence of a rectangular cylinder. 
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Fig. 13. Pattern of an electric current line source in 
the presence of a rectangular cylinder. 
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IV. CONCLUSIONS 


The GTD has been applied to calculate the radiation from a perfectly- 
conducting rectangular cylinder in the presence of a linear array of line 
sources, which may be of the electric current, magnetic current or mag- 
netic current moment type. When densely packed, these sources may be 
used to approximate the radiation from an aperture. To insure good 
accuracy in the fields calculated from the solution described here, the 
separation of source and field points from the edges of the cylinder 
and the separation of the edges from each other should be not less than 
0.7 wavelength. However, for far-zone pattern calculations, the line 
sources can be only a few tenth of a wavelength from the nearest edge. 

The use of new scalar diffraction coefficients valid in the transi- 
tion regions makes it possible to calculate continuous patterns in the 
region surrounding the cylinder away from its edges. Radiation patterns 
calculated from this solution and from an integral equation solution are 
found to be in excellent agreement for a number of stringent test cases. 
This demonstrates the utility and accuracy of the new diffraction 
coefficients and the overall accuracy of GTD as it has been applied to 
this problem. 


26 



APPENDIX I 


THE FIELD AT THE SHADOW BOUNDARY OF A 
THICK SCREEN FOR GRAZING INCIDENCE 

In this appendix we derive an expression for the field near the 
shadow boundary of a thick, perfectly conducting screen illuminated by 
a li he source at grazing incidence, as shown in Fig. 14. The solution 
near the shadow boundary in the forward direction is of interest. In 
the following development we employ Eqs. (5) through (10) in the text, 
the subscript h on the hard scalar diffraction coefficient has been 
omitted, and it is convenient to use the function 


-jkx 

(A-l) f(x) = e —_ 
vx 



Fig. 14. Shadow boundary of a thick screen for grazing 
inci dence. 
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Let the strength of the line. source be such that the incident field 
at P is 


(A-2) U i (P) = f(fc 0 ) . 

The total field at P is the sum of the inci dent field plus the field 
of the ray singly-diffracted from edge 1 and the field of the ray doubly- 
diffracted from edge 2. In the illuminated region, e > 0, 

(A-3a) U(P) = f ( £ 0 ) + f(£ i ) f(^) + HU j (2) 0(^,0, L) H^) . 

In the shadow region, e < 0, 

(A-3b ) U(P) = HU 1 (2) D(<j> 2 ,0,L) f(A 2 ) . . 

Here, 

(A-4a) <j> 2 = ir - e , 

A,* An 

(A-4b) L x = — , 

£ i + '*i 

(A-4c) U 1 (2) = fUi+h) + f(A.j) - |e J ^ cot(~) F(2kLp f(h) 
in which 

A- h 

(A-4d) L; = — ^ — 

Aj + h 

Thus for e > 0, 

(A-5a) U(P) = f(A 0 ) + f(A i ) f(A 2 ) + ^(A^h) D(<j> 2 ,0;L 2 )f(£ 2 ) 

f o J 2 it 1 

+ J af(A i )f(h) - j e cot(y-) F(2kL ’ ) Df^.O.E) f(i 2 ) , 
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and for e < 0, 


(A-5b) U(P) = hfUi+h) D(4» 2 ,0,L 2 ) f(t 2 ) + %f(£ i )f(h). 

2 "J 2tt 

• - cot(— ) F(2kL j ) 

where 

U,- + h) £ 2 
(A-6) L 2 = •— ! 

£.j + h + £ 2 

and L is a distance parameter determined by the wavefront curvature of the 
field incident on edge 2 which has been singly-diffracted from edge 1. 
Since edge 2 is in the transition region of this field, the curvature of 
this wavefront is not simply that of a cylindrical wave emanating from 
edge 1; i.e., E ^ h £ 2 /(h+£ 2 ). We will determine E by requiring U(P) to 
be continuous at the shadow boundary. 

As e -> 0, 


D(4» 2 ,0,L) f(£ 2 ) , 


(A-7) D(<t> 2 ,0,L ) = 


_ e -J it/ 4 
n^lT 


(I s -) 


iTkL2 j tt/ 4 
« Id'e 


= - vEJ sgn e 

In a similar manner, 


( A-8) D( <}> 2 ,0 ,ET) = - vC sgn .e 


Furthermore, as e -> 0 


(A— 9) 





F(2kL ! ) 

/2irk 


Substituting Eqs. (A-7), (A-8) and (A-9) into Eqs. (A-5a), (A-5b) and 
requiring U(P) to be continuous at the shadow boundary e = 0, we see that 
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(A- 10) 


e jk(Jl i +£ 1 ) 




cotCj 1 ) 


F(2kL x ) 
/2irk 


-jk(«,j+h+)t 2 ; , 




7T / 4 


2, f ( 2kL l> 

3 .'KiT 



/f 


where h + i 2 = Aj at e = 0. 
From which 


(A-ll) 

r- f h H V 

V> + V 

wi th 


(A-12) 

m = 

F(2kL;) 
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APPENDIX II 


DESCRIPTION OF THE COMPUTER PROGRAM 


A. Input Variables 

N : is the number of sources in the array which approximates the 

aperture distribution. Here N = 2M+1, where M is an integer 
(.. N is an odd integer), and M has been introduced earlier in 
section II. The DIMENSION cards at the beginning of the 
program must be dimensioned as N or larger. 

TYPE : is a reference parameter. TYPE is set equal to 1.0 when sources 

of type I (see section II) are used. TYPE is likewise set equal 

to 2.0 for type II, and is set equal to 3.0 for type III sources, 
respectively. 

AL : is the aperture width (= W of Fig. 7). 

AM(I): is the magnitude of the Ith source in the array which approximates 

a given aperture distribution. 

AP ( I ) : is the phase of the Ith line source in the array (which approxi- 

mates a given aperture distribution), in RADIANS. 

X : is the point of incidence on the 2-D box and corresponds to X 

shown in Fig. 7. 

XL : is the length of the box (corresponding to XL of Fig. 7). 

H : is the height of the box (it corresponds to h in Fig. 7). 

RI : is the incident range (corresponding to R-j in Fig. 7). 

RS : is the scattered range (corresponding to R s in Fig. 7). 

THI : is the angle of incidence (corresponding to 0^ in Fig. 7). 

XLAMDA: is the transmi tted wavelength . 

Note that the variables AL, XL, H, X, RI and RS have the same units as 

XLAMDA. 

B. Output Variables 

THS : is the angle of scattering (corresponding to e s in Fig. 7) in 

degrees . • 

ATAL : magnitude of the total field. 
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DBTAL: magnitude of the total field in dB. 

DBSAS: magnitude of the scattered field in dB. 

DBGA : incident field with phase center at Q (Fig. 7) as a function of 

THS , in dB. 

PHASE: phase of the total field in degrees. 

TH : angular variable corresponding to (Fig. 7), in degrees 

DBHA : incident field with phase center at 0 (Fig. 7) in dB, as a 

function of TH. Note that DBHA is an output variable in the 

Subroutine TEST. 

C. Instructions for Representing Aperture Field Distribution by a 
2-D Line Source Array. 

When dealing with the input variables AM(I) and AP(I) for Ith source 
in the planar array used to approximate a given aperture distribution, 
the ordering of the array elements is done as follows: 



Fig. 15. The ordering arrangement of the line sources used 
to approximate a given aperture distribution. 


Let N 8 11) where N is the number of line sources approximating a 
given aperture distribution over the aperture width- W. The ordering 
arrangement for these sources is indicated in Fig. 15. The source at the 
center of the aperture is the one for which 1 = 1. 1 = 2, 3, 4, 5 and 6 for 

sources to the right of the one designated by I = 1 (as one views the 2-D 
box from the aperture center). Similarly, I = 7,8,9,10 and 11 for 
sources to the left of the source at the center (designated by I = 1). 

Let the aperture distribution (assumed known) be represented by the 
quantity F = iFje 1 ^ over the aperture. |F| represents the magnitude of 
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the field distribution over the width W, and > represents the phase of 
the field distribution over the width W. Hypothetical plots of |F| and 
ip over the aperture are indicated below: 


IFl 



AM (I) 
AM ( 9) 

AM ( 5 ) 



AP ( 9) 

AP ( I ) 
AP( 5) 


Fig. 16. Field distribution over the aperture. 


Fig. 16 clearly indicates the amplitudes and phases of the sources 
designated by I = 1, 5 and 9. For example, the magnitude of the line 
source strength corresponding to 1 = 9 is given by AM(9), and its phase 
is given by AP (9) . Similarly, one can obtain the amplitudes and phases 
of all the other line sources. Note that the aperture is divided into 
2M segments, where M = 5. Hence, the number of sources, N = 2M+1 = 11. 
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D. Instructions for Using the 'Obliquity Factor 1 

The fields radiated by apertures are non -symmetrical on either side 
of the aperture, in most practical cases. The fields radiated by the 
2-D line source array discussed above are symmetrical on either side 
of the planar 2-D array. Thus, an obliquity factor of the type cos n e/2 
(please refer to the discussion in section II) is included for computing 
the field radiated by each source in the array. The obliquity factor is 
different for each n, where n = 0,%,1,2. A function statement FB(SX) 
computes this obliquity factor for a given value of n. Specifically, the 
statement concerning FB(SX) reads: 

FB(SX) = ABS(C0S(SX/2 .0) )**2.0 

and corresponds to an obliquity factor with n = 2.0. If any other value 
of n is desired, the appropriate value must be punched into a new card 
which replaces the previous one. Note that the value of n directly 
follows the ** symbol in the statement. 

The obliquity factor cos n e/2 is plotted as a function of e for 
different values of n (n = h, 1 and 2) in Fig. 17. When the pattern of 
an isotropic source is multiplied by cos 11 e/2, it is evident from the 
resultant pattern that the obliquity factor serves to control the level 
of the radiation pattern primarily in the range tt/3< e< 5tt/3 . The case 
n = 0 corresponds to the isotropic case. 

E. Instructions for Computing the Incident Field 

Two incident field patterns are computed, one .is for a phase reference 
at the center of the aperture, and the other is for a phase reference at 
Q (see Fig. 7). In the former case, the radiation pattern in dB is desig- 
nated by DBHA, and is obtained as a function of e (or TH as defined in 
the computer program). In the latter case, the radiation pattern in dB 
is designated by DBGA and is obtained as a function of e s (or THS as 
defined in the computer program) . DBGA is computed at a distance equal 
to R s from Q. DBHA has been programmed for a range of R-j + R s from 0 
(center of the aperture as shown in Fig. 7); however, if the user wishes 
to change the present range for DBHA , only one card in the program deck 
needs modification. A subroutine designated TEST computes DBHA at a 
range of Ri + R s from 0; the call statement for this subroutine is 

CALL TEST (N, AL, RS+RI , A, TYPE) 

If a different value of the range is desired, one must replace RS+RI in 
the call statement above by a number which equals the desired value for 
the range. Note that the new range should have the same unit as those of 
X (corresponding to XL AM DA in the computer program). 
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Fig. 17. Patterns of the obliquity factor for ditferent 
val ues of n. 


F . Sample Programs 

In this section, we present a sample case which serves to illustrate 
the use of our computer program. The example selected involves an array: 
of three magnetic line sources of unit strength which illuminate a rec- 
tangular cylinder, as in Fig. 18. We utilize the computer program for 
calculating the incident field of the array, the field scattered by the 
2-D box (rectangular cylinder) and the total field (incident + scattered) 
surrounding the 2-D box. 


35 









Fig. 18. An array of three line sources in the presence 
of a square cylinder. 


The ordering of the array elements is shown in Fig. 18, where the source 
at the center is labeled source #1. Note that source #1 corresponds to 
1=1, and sources #2 and #3 correspond to I = 2, and I = 3, respectively. 
For this particular problem, AM(1), AM(2) and AM(3) are each equal to 
1.0, and AP ( 1 ) , AP (2 ) and AP(3) are each equal to 0.0, because the line 
sources are of unit strength and zero phase. An obliquity factor corres- 
ponding to n = 2 (i .e. , obliquity factor = cos 2 e/2} has been incorporated 
into the program for the incident field pattern, and the incident field 
pattern corresponding to DBHA (phase reference at source #1) is plotted 
in Fig. 19. Also included in Fig. 19 is the incident field pattern with- 
out the obliquity factor (n = 0 case) for the sake of comparison. The 
pattern of the scattered field designated by DBSAS, and computed for 
values of e s (or THS) which lie in the range -180° < 0s i 180°, is plotted 
in Fig. 20. The scattered field obtained by our method is compared 
against that obtained from a numerical solution to the integral equation 
for this problem given by J. H. Richmond; these results agree perfectly. 
Finally, the total field (incident plus scattered, each being phase 
referenced at Q) is also obtained, and is designated by DBTAL. DBTAL is 
computed as a function of e s (THS in the program) and the results are 
indicated in Fig. 21 by a dashed curve. The solid curve is added for the 
sake of comparison; it corresponds to the total field when the incident 
field has no obliquity factor in it. 
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WITHOUT THE 
OBLIQUITY FACTOR 



180 ° 


Fig. 19. Patterns of an array of three magnetic line sources 
of equal strength, with and without the obliquity 
factor. 
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Fig. 20. Pattern of the field scattered by a square cylinder 
which is illuminated by an array of three magnetic 
line sources. 
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FOR ILLUMINATION 
WITHOUT OBLIOU IT Y 
FACTOR 



180 ° 


Fig. 21. Pattern of an array of three magnetic line sources 
in the presence of a square cylinder. 
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